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Abstract 



The matrix model with a Bethe-tree embedding space (coinciding at large 
iV with the Kazakov-Migdal "induced QCD" model |]J) is investigated. We 
further elaborate the Riemann-Hilbert approach of Ref. assuming certain 
holomorphic properties of the solution. The critical scaling (an edge sin- 
gularity of the density) is found to be j str = — ^arcosD, for |D| < 1, and 
Istr = — ^arcos for D > 1. Explicit solutions are constructed at D — ~ 
and D = oo. 



1 Introduction 

Let us start with the matrix model describing surfaces embedded into an 
arbitrary graph Q defined by its incidence matrix G xy {G xy — 1, if there is a 
link between vertices x and y, and 0, otherwise). The standard action is 

S = -N G xy ti [$(£)<%)] + N tr V($(x)) (1) 

xy x 

where <&(x) is an iV x N hermitian matrix attached to an x'th vertex; the 
potential, V(ip), is an arbitrary polynomial. The partition function is defined 
as the integral over all field configurations 

z = I n d ^ x ) e ~ s ( 2 ) 

J x&g 

At each vertex, can be decomposed into diagonal, (p(x), and angular, 
S(x) E U(N), parts: 

= S + {x)p{x)S{x) (3) 
It is convenient to introduce the on-link gauge variables 



Q xy = S + (x)S(y) (4) 

obeying the constraint that their product along any loop equals unity. 

If the graph Q is a tree, all the gauge variables are independent and can 
be integrated out by the Itzykson-Zuber formula 0: 

J A(0)A(^) 

In eq. (^), (j) and ip without loss of generality are real and diagonal; 
A(0) = rii<i(0j — 4>j) is the Van-der-Monde determinant. 

As a result, one finds a model in which a role of dynamical variables is 
played by iV eigenvalues of In the iV — ► oo limit, the saddle-point 

approximation is exact and, in this sense, the model is soluble. 

If the graph Q has loops, the model becomes very complicated, and so far 
no reliable results have been obtained in this case. Technically, mean field 
does not work; physically, models of this kind describe the interaction of 2d 
gravity with matter having the central charge c > 1. 
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In the simplest case of the c = 1 model compactified on a ring, the angu- 
lar degrees of freedom describe vortices, i.e., configurations with holes in the 
string world sheet wrapped along the ring. By dropping the constraint that 
the product of the angular variables along the ring equals unity, one discards 
the vortices [3J. Hence, the singlet (vortex free) sector of the model is phys- 
ically sensible. By dropping the constraints in the matrix model embedded 
in the regular D-dimensional lattice, Kazakov and Migdal Q obtained the 
model (sometimes called "induced QCD" ||) generalizing the c = 1 vortex- 
free matrix model. It can be easily seen that, in the N — > oo limit, it is 
equivalent to the matrix model with a Bethe-tree (BT) target space ||. A 
Bethe tree is, by definition, an infinite tree having the same even coordina- 
tion number, 2D, at each vertex. It is the infinite simply connected covering 
of the regular D- dimensional lattice. Therefore, the BT matrix model can 
be regarded as the mean field approximation for lattice scalar field theory. 



2 The matrix model on a Bethe lattice 

In this paper we consider the matrix model with a BT target space. In the 
large N limit, the ground state should be homogeneous. The free energy per 
lattice site is given by the matrix integral 



F = log / d N X e~ Ntvy W [l(X) 
where I(X) obeys the equation 



I(X) 



d N 2 y e Ntr(XY-V(Y)) 



I(Y) 



2D 



2D-1 



In terms of eigenvalues of X we find 



r N N 

F = log / J] dx n A 2 {x) e~ N ^i v ^[l(X) 

^ n=l 

The saddle-point equation for this integral reads 



2D 



(6) 



(7) 



(8) 



E 



./V ■ i • Ob i X n 



— V'(x i )+2Dw i = 



(9) 



2 



where 

Let us assume that, in the N — > oo limit, the eigenvalues distribution 

1 N 

ra=l 

has a finite connected support. Then we can introduce the function 

m=ldyML (12) 
J x — y 

where the integral goes over the support of p(y) (Supp p{y)). f(x) is holo- 
morphic at the infinity f(x) = - + O(^) and has one cut on the real axis 
where 

Im/(x) = iip(x) (13) 
In terms of f(x), Eq. @ takes the form 

2Re f(x) - V'(x) + 2Dw(x) = (14) 
which is valid, strictly speaking, only on Supp p(x). 

w(x) = T x ^)^W losI{x) (15) 

is the large iV limit of the function (|T0|). 

To close the system of equations, we should calculate w(x) as a functional 
of f(x), i.e., to solve Eq. (|7|). However, we can avoid doing these complicated 
calculations. As was shown in Ref. (see also 0), in the iV — > oo limit, 
the following equation holds 

l = (z- w(x))W(z, x) - j dyp{y) W ^' x) ~ W{ - Z ' v) (16) 

J x — y 

where W(z, x) is the N —>■ oo limit of the mean value of diagonal elements 
of the resolvent matrix 



W k (z) 



I(X) lz 



1 



3 



I(X) (17) 

kk 



The large z expansion of this function has the form 

• , 1 w(x) ,„ . 

W (z,x ) = - + — V Z + ... (18) 

z z z 

As long as the ground state is considered, the following homogeneity 
condition holds 

' dx p{x)W{z,x)= f{z) (19) 

Eq. fll6|) is valid only on Supp p(x). However, all functions can be ana- 
lytically continued into the whole complex plane. Following Ref. 0, let us 
introduce the function 



p(y)W(z,y)^ 1 _ e _ F f dN , xdN2 



F(z, x) = l- / dw— — = 1 - lim e~ F d N X d N Y 

' X-y N^oo 1 



1 

— tr 



1 1 



c iVtr(XV -V'(X) -v(y)) - - 2i> 1 



v v AX)I(Y) (20) 

_z — X x — Y \ 1 J 

where the integral goes over Supp p(x). By construction, F(z, x) is symmetric 
F(z, x) = F(x, z). 

At z sufficiently large, F(z, x) as a function of x is holomorphic every- 
where except the support, where its imaginary part equals 

Im F(z, x) = -irp(x)W(z, x) (21) 
The real part can be determined from Eq. ( |I6|) 

Re F(z, x) = (z- w(x) - Re f(x))W(z, x) (22) 

As F(z, x) is real on the real axis out of the cut, we find from the Riemann- 
Schwartz principle that 

F(z,x) = F(z,x) (23) 

and, hence, the ratio of limiting values of F(z, x) above and below the cuts 
is independent of W(z, x) 

F(z,x + i0) z — w(x) — Re f(x) — ilm f(x) , . 

F(z,x — iO) z — w(x) — Re f(x) + ilm f(x) 
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At the infinity we have by construction 



F(z,x) 



1 - 



/(*) 



X y x z 

The function having all needed properties is given by the formula 



F(z,x) 



exp 



dy 



log[z - w(y) - f{y)] 



(25) 



(26) 



Jc 27rz y — x 

where the contour C encircles the cut of f(y) in the positive direction leaving 
aside all other singularities (whatever they are), that is, F(z,x) is defined 
perturbatively by its expansion in inverse powers of z and x. 

For Eq. (p6|) to make sense, the function u(x) — w(x) + f(x) has to be 
holomorphic. However, w(x) depends through the saddle-point equation ( P|) 
on Re f(x). Hence, we should be able to introduce two holomorphic functions 
fi(x) and f2(x) such that, on Supp p(x), fi(x) = Ref(x) and f2(x) = 
lmf(x). Off the support, both fi(x) and f2(%) are holomorphic (i.e., obey 
the Cauchy-Riemann equations separately); on it Im/i(x) = Im/ 2 (^) = 0. 
Everywhere in the complex plane f(x) = fi(x) + if^x), and we define 



<y) = ^V'(y) + ^f 1 (y) + iMy) 



(27) 



Let Supp p(x) = (a, b), then we can rewrite Eq. (|26D as 



F(z, x) = exp — / 

J a 



dy 



2m y — x 



log 



z - u(y) 



u(y) 



exp 



b du 



2ni u 



log 



x - y[u) 



z 



y{u) 



(28) 



where y(x) is the inverse function: u(y(x)) = x. For F(z, x) to be symmetric, 
F(z, x) = F(x, z), the following remarkable equation has to hold 



y[x) 



u(x) 



(29) 



or, equivalently, 



u(u(x)) = X 



(30) 
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This equation has previously been obtained by Luochen (unpublished) and 
rederived by Matytsin within a quite different approach ||. 

In general, the inverse function y(x) is not unique. However, we are 
interesting in critical behavior consisting in the divergence of a perturbative 
expansion in a non-gaussian part of the potential V(x). In the gaussian case, 
V (x) 



^-x 2 , the solution was found by D. Gross 



/(*) 



a / a 2 x 2 

—x + i\ a 

2 V 4 



with 2 possible values of a 

1 



m\D - 1) ± DJm 4 - 4(2D - 1 



(31) 



(32) 



1 2D - 1 

The condition that a solution has to describe perturbations around this semi- 
circular saddle-point fixes a branch of y(x) unambiguously. 

From the expansion (P5D, we find Master Field Equation in the form 

/(*)=/ ^-\og[z-w{y)-f{y)\ (33) 



which is satisfied automatically, if Eq. (|30|) holds. 
We can as well obtain the equation 



/(*) 



—^r log z — w(y) — f(y) 



then, obviously, 
dy 



lc 2iri 
and we conclude that 



{f2(y)) : 



(34) 



(35) 



f 2 (x) = yj$(x) (36) 

where $(x) is holomorphic inside the contour C. Critical behavior takes place 
when (i) a branching point or (ii) a zero of $(x) approach the support of the 
eigenvalues density (its edges correspond to zeros of <&(x)), in other words, 
when two singularities pinch the contour. In the first (i) case, the order 
of the branching point increases, which is characteristic of c < 1 models. 
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The second (ii) possibility corresponds to the shrinking of a handle of the 
Riemann surface of f(x), as in the c = 1 matrix model. 

3 Critical scaling 

In matrix models, critical behavior is determined by a scaling near an edge 
of the eigenvalues density exactly at a critical point, when 

p(x) oc (-x) 1+7 x < (37) 

where 7 = — j s t r > is a critical exponent (for simplicity, we have assumed 
that the branching point of f(x) is at x = 0). In the small vicinity, \x\ <C 1, 
we can expand all functions in powers of x. For example, (if < 7 < 1) 

u(x) ~ ax + — — — cos 7T7 + % sin 717) e - ™ 1 x l+1 + . . . (38) 

and 

u{x) Rj ax + — — — cos 7T7 — z sin 717) e _i7r7 x 1+7 + . . . (39) 
Expanding Eq. flSTf), we find 



x = m(-u(x)) ~ a 2 x + abe l7r Y x i+i cos 717 + i sin 7r7 + {^jp cos 717 — z sin 7^ a 7 

(40) 

Hence, a 2 = 1. If we take a = 1, a 7 = 1, the only solution is 7 = |, which 
is the well-known one-matrix-model exponent. Indeed, this solution should 
always be possible. For example, at D — |, it corresponds to the disordered 
phase of the Ising model. 

The second possibility is a = —1, which can be interpreted as a quadratic 



x 



top of the effective potential, V e ff = j^V'(x) + ^rfi e9 (x), where f[ 
is a regular at x = part of fi(x). Therefore, one should expect critical 
behavior in this case. In principle, we can take 2 branches of the function 
a 7 . In the first case, a 7 = e 171 " 7 , we find the equation 

cos7T7 = D (41) 

which makes sense for \D\ < 1 and fits 3 exactly soluble cases: D = 0, | and 
1. 
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The second possibility is a 1 = e 11X1 . In this simple algebra gives 

the equation 

COS7T7 = — — ^ — (42) 
which formally makes sense for D > 1 and D < |. 

4 Explicit solution at D = | and oo 

In the previous section we have determined possible types of singularities of 
p(x). However, it would be instructive to have a non-trivial global solution 
to the problem. 

Let us suppose that there exist holomorphic coordinates A in which the 
function f(x) becomes an integral function: 

/(:r(A)) ee (^(A) = ]>> fc A fc (43) 
fc>i 

and the infinity of the x-plane is mapped to the center of the A-plane (i.e., 
x ~ j-, when A — > 0). It means that, geometrically, the Riemann surface of 
f(x) is a sphere with one puncture. 

There should exist an involution transformation of the A-plane, 7r(A), 
preserving a curve of real x's, 

7r(7r(A)) = A x{n(X)) = x(X), (44) 

such that, on the curve (and only on it), the action of 7r coincides with the 
complex conjugation: 



<p(\) = (p(n(X)) (45) 
Then, as f(x) = f(x), we can identify 

^i(A) = / 1 (x(A)) = |b(A)+ ¥ »(7r(A))] 

(46) 

W ! (A) = / 2 (x(A)) = ib(A)- ¥ »(7r(A))] 
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The functions u(x) and u(x) define an invertible map of the A-plane onto the 
/i-plane and back, correspondingly, determined by the equations 



x(jm) = ±V'(x(\)) + ^V(A) - ^(tt(A)) 

x(X) = ±V'(zQ*)) - ^tp(n) + ^M/,)) 
The simplest realization of this scheme, 



(47) 



x(X) = - + a\ tp(\) = A, (48) 

A 

corresponds to the semi-circle solution fl3~l|) . Here, 

tt(A) = -±- (49) 

CLA 

For the gaussian potential, V'(x) = m 2 x, equations (|47|) take the form 

i + aA = 2^(m 2 a + 2D - 1)± + ^(m 2 a - 1)// 
from which we find 



(50) 



/i = — -A 51 

m^a — 1 



together with the self-consistency equation for a: 



(m 2 a + 2D- l){m 2 a- 1) = (2Da) 2 (52) 



2 



which gives the solution (f?2j). 

Now, let us consider the non-gaussian potential ^V'(x) = m 2 x + gx 
and try the ansatz 

x(X) = \ + c + aX + bX 2 (53) 
A 

From the power counting, we find the form of tp(X): 

(p(X) = A + ip 2 X 2 + ip 3 X 3 + tp 4 X 4 (54) 
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In this case, 7r(A) obeys the quadratic equation 




b\iT 2 (X) + (aX + 6A 2 )tt(A) -1 = (55) 
which has two solutions 

7r ± (A) = ~[-^-A±^ + (~ + Ar (56) 

We have to use 7r + (A) for the conjugation, because it is an involution: 7r + (7r + (A)) 
A. Its asymptotic behavior is 

f -?=, when A — > 
tt + A = fA' 57 
I 6A2' w hen A — > oo 

It is natural to suppose that zero and the infinity are fixed points of the 
function /i(A) (/i(0) = and /i(oo) = oo) as it was in the gaussian case. We 
see then that a number of terms on the right hand side of Eq. (^) have to 
cancel identically, which gives us the asymptotics of /i(A): 



f\\ — I Vy^' when A — > , , 

M } " 1 v%A 2 , when A - oo (58) 

together with the answer for tp(X): 

ip(X) = A + [m 2 b + g{a 2 + 2bc)]X 2 + 2abgX 3 + b 2 gX A (59) 

If D ^ 2> we find contradictions and the equations have no solution. One 
could have expected it from the very beginning, because the ansatz ( |53| ) is 
well known in the orthogonal polynomials approach to the two-matrix model 
and gives -y str = -\ 0. 

At this point, Eq. fl47j) take the form 



(60) 



j^ + c + a^i + bfi 2 = b 2 gX 4 + 2abgX 3 + (m 2 b + g(a 2 + 2bc))X 2 + 

((m 2 + 2cg)a + bg)X + m 2 c + (2a + c 2 )g+ 

{b 2 gX 3 + abgX 2 + (m 2 + 2cg)bX + 2&a + (m 2 + 2 C£ /)a - 1}tt + (A) 



\ + c + aX + bX 2 = \ + [m 2 + 2gc)- + (2a + c 2 )g + m 2 c + (61) 
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Let us assume that the function /i(A) obeys a quadratic equation, then 
there is no ambiguity left and the ansatz having correct asymptotic behavior 
is of the form 

fji(X) = -Jbg\n_(\) = Jbg\7t + (X) + M{aX + b\ 2 ) (62) 

It can be easily checked that this function does satisfy Eqs. ( |60| , |6TD provided 
the parameters obey the following three constraints 



a 



Mm 2 + 2gc) 



gc 2 + (m 2 - l)c + 2ag = (63) 
2bg - + (m 2 + 2gc)a -1 = 

which completes the solution in the D — A case. 

If D = oo, Eq. (|47| ) takes a similar form and the same ansatz fl53|) should 
work. In this case, it is convenient to introduce the function A(/i) which 
coincides with Eq. (|62|): 



A(/x) = -jbgiw.bi) (64) 
and, in complete analogy with the previous case, we find 

<p(\) = A - \m 2 b + <?(a 2 + 26c)]A 2 - 2abg\ 3 - b 2 gX 4 (65) 
where parameters obey the equations 

a = ^(m 2 + 2gc) 

gc 2 + (m 2 - l)c + 2a# = (66) 
2bg- (m 2 + 2gc)a + 1 = 

The only difference between Eq. (|63"D and Eq. (j66|) is the sign before 1 in the 
last equation. 

At the critical point we have 



(1 + sA) 3 

x(X) = V J (67) 



and a simple algebra gives 



s2 = | ^ whenD = | 
1 — - when D = oc 



lo- 



ll 



The critical coupling, g c , is to be determined from the equation 



Vs9~c = 1 (69) 

Hence, at D = oo, g c is complex. 

5 Discussion 

The main result of this paper, Eq. (PI) , looks very natural. After the iden- 
tification D = |, it coincides with the Kostov's solution of the O(n) matrix 
model |TI|. It is a standard matrix-model singularity. The second branch 
( (42"|) corresponds to a perturbation expansion around a local maximum of the 
free energy Q (only this saddle-point of the gaussian model is non-singular 
at D = oo). Moreover, we have found that the critical coupling in this case 
is complex (at least, at D = oo). This means, presumably, that, at D > 1, 
our model is in the usual sense non-critical as all solved so far matrix models. 

6 Acknowledgments 

I would like to thank J.Ambj0rn, V.Kazakov, Yu.Makeenko and, especially, 
A.A.Migdal for the discussions. Financial support from the EEC grant CS1- 
D430-C is gratefully acknowledged. 

References 

[1] V.A.Kazakov and A.A.Migdal, Nucl.Phys. B397, 214 (1993). 

[2] A.A.Migdal, Mod.Phys.Lett. A8, 359 (1993). 

[3] C.Itzykson and J.-B.Zuber, J.Math.Phys. 21, 411 (1980). 

[4] D. Gross and I.Klebanov, Nucl. Phys. B344, 475 and B354, 459 (1990); 
D.V.Boulatov and V.A.Kazakov, Int. J. Mod. Phys. A8, 809 (1993). 

[5] A.A.Migdal, Mod.Phys.Lett. A8, 153, 245 (1993); 

Yu.Makeenko, Mod.Phys.Lett. A8, 209 (1993), preprint ITEP-YM-4-93 



12 



(1993); 

V.A.Kazakov, Nucl.Phys. (Proc. Suppl.) 30B, 149 (1993). 
[6] D.V.Boulatov, Mod.Phys.Lett. A8, 557 (1993); 
[7] D. Gross, Phys.Lett. 293B, 181 (1992). 

[8] A.Matytsin, On the large N limit of the Itzykson-Zuber integral, Prince- 
ton preprint (1993), hepth/9306077. 

[9] J.-M.Daul, V.A.Kazakov and I.K.Kostov, Rational theories of 2d gravity 
from the two-matrix model, preprint CERN-TH. 6834/93 (1993); 
M.Douglas, in Proceedings of 1990 Cargese Workshop on 2d Gravity. 

[10] I.K.Kostov, Mod.Phys.Lett. A4, 217 (1989). 



13 



